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I. INTRODUCTION 
A. PURPOSE AND OBJECTIVE 

The purpose of this thesis was to investigate a method 
of moment comparisons to determine statistical descriptions 
of aircraft structural fatigue characteristics through 
appropriate selection of probability distribution functions. 
Once the population distribution and it's parameters are 
determined, they can be used to model fatigue parameters 
such as applied loads that are damaging, static and fatigue 
lives of material characterization samples, full scale 
fatigue tests, air gust load intensity, sink speeds for 
carrier landings, etc. 

Realities are such that in many of these instances the 
sample size may be rather small, statistically speaking. 
The full scale fatigue test is a sample size of one because 
of cost, and yet some estimates about other aircraft in the 
population must be made from it. This represents’ the 
extreme, but there are many instances in helicopter flight 
Substantiation tests where the number of load variations in 
a maneuver will be on the order of 20. For this reason, a 
second objective of this thesis was be to study sample size 


influences on the methods of characterization. 


B. BACKGROUND 
Fatigue lives are based on two factors: stresses 
endured by the aircraft, and the material fatigue 
properties. Both of these factors must be studied through 
statistical analysis due to random variations in their 
values. 
1. Statistical Nature of Flight Loads 
In aircraft fatigue work, fatigue damaging stresses 
result from loads encountered by the aircraft being studied. 
In the case of helicopters, these loads are measured during 
the final design phase through substantiation flights. 
Critical components and structural members are fitted with 
strain gages, and a time-load history 1s _ recorded. 
Typically 40 to 60 parameters are recorded in each of a 
variety of flight regimes which reflect the expected or 
potential usage and flight envelope of the aircraft. An 
effort is made to have the same gross weight, center of 
gravity, airspeed, rotor speed, load factor, density 
altitude etc., for each flight to minimize variabiliGyaeee 
the loads produced. However, items like pilot technique, 
gusty air conditions, and instrument accuracy introduce 
considerable variations. Because of the high cost of 
Substantiation flights, it is not feasible to gather large 
amounts of data by flying hundreds of flights to 


statistically establish the "true" loads encountered in each 


flight regime. One method used to ensure conservatism is to 
select the highest load encountered in a given regime, and 
assume it occurs 100% of the time the aircraft flies that 
regime. Another method used is to select the 95th 
percentile load of those recorded in a given regime and 
assume it occurs 100% of the time the aircraft flies that 
regime. Both of these methods are likely to be 
conservative, but they are also subject to "high envelope 
growth". This occurs when subsequent flights produce higher 
loads than the established highest load (or a higher 95th 
percentile load). [Ref. 1] When this occurs it reduces the 
validity of the fatigue life predictions which were based on 
earlier substantiation flights. It would be useful to 
determine the population distribution of the load samples 
from each flight regime. Statistical models could then be 
constructed to model the loads for situations not covered by 
the substantiation flights, or for situations that were not 
represented accurately during substantiation due to the 
uncontrollable factors mentioned earlier. 

Every Navy aircraft is monitored fatigue-wise. Many 
are monitored by making load measurements using a variety of 
instruments. Measured data is recorded in the fleet and 
sent at regular intervals to designated ground facilities 
for processing. As a result of missing data, the gaps are 


filled via statistical representations of loads data from 


all the other aircratt (of) Chaeeeeye ce This statistical 
application to the fatigue calculation is dependent upon the 
Size of the fleet of that particular aircraft, which is 
sometimes small. 
2. Statistical Nature of Material Behavior 

Even if the stresses or loads to be encountered were 
known exactly, fatigue life still could not be calculated 
exactly due to the inherent variations in the strength of 
the structures or components themselves. Manufacturing 
processes like stamping, rolling, grinding, machining or 
heat-treating can cause surface irregularities or residual 
stresses. Assembly processes such as rivet placement and 
bolt torque can produce variations in stress concentrations. 
Even nondominant, microscopic, interstitial impurities or 
dislocations have a strong effect on fatigue life. 
Depending on the loads involved, these small irregularities 
can cause a large scatter in fatigue lives for material 
samples, subassemblies, and full scale tests, even when 
subject to identical loading. "Coupon" testing is 


relatively inexpensive, so sizable amounts of data can be 


generated, which can. be Statistically analyzed to 
characterize the variability. If the data can be fit to a 
distribution and modeled, various "“safe-lives" can be 


developed to ensure an item will not fail before a certain 


period of time. These lives are assigned confidence levels, 


which indicate that a specified percentage of the items 
would statistically survive the calculated safe-life a 
specified percentage of the time. [Ref. 2] 

Having motivated the broad need of statistical 
characterization in fatigue analysis of aircraft, the moment 


method fundamentals will now be developed. 


II. DISTRIBUTION MOMENTS METHOD 

A. MOMENTS DEFINED 

The method of distribution identification developed in 
this chapter is based on a group of statistics known as 
moments of the probability distribution function. The first 
moment of a distribution is more commonly referred to as the 
mean [l., or the expected value of the random variabaemeas 
The expected value is the centroid of the distribution. 


M, = ELX]= | xf (x)dx ee 


—o 


BG is the expected value, or expectance. For a continuous 
distribution, f(x) is the distribution's probability density 
function, and x is the value of the random variable. The 
mean is the first moment about zero. Higher order moments 
can be calculated about any point in the distribution, but 
are most commonly taken about zero (the raw moment), or 
about the mean (the central moment). The second central 
moment up’ is known as the variance, and is a measure of the 
distribution's dispersion. It is akin to the momentwies 
inertia in dynamics or structures. The square-root Of eis 
variance is the more commonly used standard deviation, 
denoted by the symbol o. Equation 2.2 is the expression for 


the second moment about the mean. 


+ oO 


we = E(x —,) |= fe-m,) fadde 2.2 
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Equation 2.3 is the general form for the rth order, central 


moment and is defined as expected: 


+00 


wo? =E|(X- 4.) |= f(x- 4) Sa Bee 
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It is important to note that the superscript on wm, in 
parentheses is not an exponent. 
The expectance term for the variance is expanded as 


shown below: 


pi? =E|(X-y,) |= £[X?]-2H,2[X]+(u, 2.4 


where the expected value of a constant is that constant, 


El u,|=p¢, and £[2}=2. 
The third central moment is called the skewness and 
according to Equation 2.3, it specializes to Equation 2.5. 


Equation 2.6 is the algebraic expansion of the third moment. 
fie =E\(X- p,)'| ols 


po = E) X°|-3y,£|X?]+3u2E[X]- up 2.6 


This will be used in derivations to _ follow. ee 6 
distribution has a non zero third moment it is not symmetric 
about it's mean. If it is negative it possesses a left 
"tail" and is said to be skewed to the left, or negatively 
skewed. Similarly, if the moment is positive, the tail is 
to the right and the distribution is positively skewed. 

The fourth, and final, moment used in this thesis is 
called the kurtosis. This is a measure of the peakedness of 
the distribution. If the tails are long, the kurtosis is 
greater than if the tails are short. The expression for 
kurtosis is given in Equation 2.7, the algebraic expansion 


is provided in Equation 2.8. [Ref. 2 and Ref. 3] 
pi = E\(X- 4.) 2a 


us? = E[X*|-4y,£[X?|+6u7£| X?|-4 EL X]+ us 2.8 


When comparing measures of the mean, standard deviation, 
skewness and kurtosis, the above quantities are 


standardized, or non-dimensionalized as follows: 


(3) u® E|x- ey 


_ —~ Hy : 
= VB, = (ui) (z[x- yy ) 


(4) 


(4) x- ) 
= B, - —3=52— = eee | = ral —3 2 poe 


; (er (z[x- poy ) 


Notation and terminology vary widely depending on the 


Mmemerence used. Here, y is the coefficient of variation, jy, 


is the measure of skewness, and Y2 is the measure of 


kurtosis. Veneer ana sy.) can, be referred to as the 


third and forth cumulants, or shape factors as well. [Ref. 4 





and Ref. 5] The -3 is a centralizing term used when 
comparing measures of kurtosis. It comes from the fact that 
the kurtosis for a Normal distribution is 3. Some 
references use this correction, others do not. ine ee has 
4 
thesis the correction was not used. Hence, = bh = ae 


om 


The underlying approach of the distribution moments 
method is that if moment values of sample data from an 
unknown population match the moments of a _- standard 
distribution, it would be the appropriate distribution with 
which to characterize the data. This is done by plotting 
the moments calculated from the data over a template of 
moment functions calculated for various standard 


amstributions. If the sample moment values fall on, or 


near, a known distribution curve, it could indicate that the 
sample's population distribution corresponds. The 
distributions treated in this moment method are the Normal, 
Lognormal, Exponential, and Weibull. Each of these 
distributions has a unique probability density function that 
can be integrated and used to derive a function for each of 
the moments defined above. Once the moment functions are 
derived, they can be plotted versus one another, and a 
template of curves can be formed. 
B. DERIVATION OF MOMENT FUNCTIONS 
1. Normal Distribution 

Equation 2.12 is the probability density function 

(pdf) for the Normal distribution. 
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The x subscript has been dropped from the mean and standard 
deviation without creating any ambiguity. The pdf is 
substituted into Equation 2.1 and a general expression for 
the moments is developed. 

(7-H) |e 


He =E|(x- ny |= ar EXP] —— 1 2.18 
J om 20 
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To simplify this integral a variable change is needed to the 








~— 
standard Normal Zz. fet om or re then Equation 
Oo 
meee becomes: 
pu =0'f zt exp|—3 2" |e 2.14 
“N20 2 


For odd integer values of r the integral is zero, but for 


even integer values of r, Equation 2.14 can be written as: 





This iS a variation of Euler's Second integral where: 


2 "T(s)= [= exp) 32d Zao 
0 


. I , , 
EHetting s=s(r+l), and substituting into Equation 2.15 


yr9elds: (Ref. 3] 


oy”? r{F(r+1) 
(r) 
. a 


Since odd integer values of r result in pw”=0, the 


standardized measure for skewness, Equation 2.10 becomes: 


JEN 


(3) (3) 


v,=/B,=+ a. Os 


x 
3 vere 
a. 





= 0 2s 


which indicates a symmetrical distribution. Using 
properties of the gamma function, the values of yw” and p™ 


ale. 


= =6 2.19 
‘ Un 
o*2'1 5 
cs 2) =36' 2220 
TT 


Substituting these values into Equation 2.11 for the 


standardized measure of kurtosis yields: 


(4) (4) 
BS LM, 30 
Y> = DB; = A = 5 = A 
6, (wy 





Recall, the correction factor of -3 is not being used. 
2. Lognormal Distribution 
Using the random variable Z, where x=Inz is normal 
and the constants £ and 6 as parameters, where w.=€ and 
o,=0, the Lognormal probability density function for Z 


becomes: 
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] 2 
f.2)= oom, ex] linz- air for z> 0 ee 


=Q0Q, forz<0 


Substituting this equation into Equation 2.1 yields a 


general form for the expected value of the rth moment of Z. 





EZ [aes] Sil E Ja Dees 


By changing the variable to y=Inz-€; then z=e’e’, and 














dz=ee’dy, and -o<y<o [Ref. 3]. Substituting and 
fap lLitying: 
Neves) aman =) 
E|\Z |= e In(e’e® e’erdy 2.24 
2"]= [enn] Se linle’e!) 2) fret 
VALE 
*s aa =2 2.25 
3 5? 5 
ea 53 ode Das 
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At this point 6‘r*? is added and subtracted to the 
exponential to make (25? + y?] a perfect square. 


E\Z’|=e"e 





—] 
Fear om Soa (y?-28 n+ 8'?) by 2.29 
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Finally: 








l r | —| 2 
E{zZ’])= re —~5’r) ld 2.30 
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The integral factor is the same as the cumulative density 


function for the Normal distribution, which is normalized to 


one. Therefore, the general expression for the expectamee 
1S: 
] 
BIZ" |= exp|r 458°" 2.30 
and, 
] 
p= E[Z]= exp] s+38"| 2.32 


The coefficient of variation from Equat rong eee 


V2 


y =(exps?-1) 2.32 


However, unlike the Normal distribution where a general form 


(r) 


Oru was derived, here each moment must be expanded 


algebraically and substituted into Equations 2.10 and 2.11. 


For example, to calculate the measure of skewness (y,), wy” 
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(3) 


and wu must first be calculated. The necessary 


substitutions with some examples of the algebra are shown 


below. Equation 2.4, wu? = EX?) 2pE[ X]+(u)? becomes: 


i = exp(2€+257]-2exp[26+57]+exp[26+5°|= exp| 2 +25°]-exp|2é + 5°] 
2253 


fee a Similar substitution, Equation 2.6, for the third 
moment py) = £|X*|-3y,£[X°]+3u2£[X]-y}, becomes: 


>) al yee | Pee es) 252). ae rei 2 | eee 3 2 
fv met? leeca Bee? jexp|2 § 26 3exp|2 ¢ 5 Jexp) § Y eiiee 5° | 


=ex9]3 4 8 |-3ex) 3 +2 8 |+3exp|3 E+ 8 |-exo]3 a e| 


2.34 
=ex9]3 += 8 |-3ex| E+ 5 |. dex é+= a 
Note that £[X°*|# py’. 
The measure of skewness, Equation 2.10 simplifies to: 
exp| 30° |—3exp|d° |+ 2 
ee 3 -19)0" | +2 a 


(exp[6?] ~ 1) 


iiesmeasure of kurtosis y,, can be derived once the 


is expanded and simplified. The algebra 


expression for yu 
is too extensive to show. Only the simplified function is 


Given. 


i) 


Y= exp| 46° |+ 2 exp|36°]+3exp|26*|-3 2.36 


The functions in Equations 2.32a, 2.35, and 2.36 Show jem 
the coefficient of variation and the measures of skewness 
and kurtosis depend only on the value of 6, the standard 
deviation of the population. 
3. Exponential Distribution 
The probability density function for the exponential 


adistribution is: 


f.(x) = Oexp(—Ox) , ifx20 2.37 


Substituting into Equation 2.la yields the geénemau 
expression for the expectance: 


E[X’ |= [x (6 exp(—@x))dx 2236 


, ; V 
Performing the variable substitution: Va and dv = @dx, 


The expectance becomes: 


ela] |(5) enw =a very 2.39 


This integral has the same form as Euler's Second Integral, 


called the gamma function [Ref. 3]. 


Mays [rear 2.40 
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Substituting v for t, and letting r+l=z, the simplified 


expression for the expectance is: 


T(r+1) 
9° 





E[X’]= 


Using Equation 2.41 and properties of the gamma function, 


(2) 


expressions for yw, yp”, yp, 


(4) 


and pp’ are shown below. 


(2) 


p22) 1 
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ot 2.43 
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Constructing the numerators and denominators for Equations 


O- 




















2.10 and 2.11, using the results of Equations 2.42 through 
2.45, produces a coefficient of variation of y=l, a measure 
of skewness of y,=2, and a measure of kurtosis of y,=9. 


These will plot as points on the template. 


di 


4. Weibull Distribution 
The parameters of the Weibull distribution are, 8B 
which determines the shape, 1/a which plays the same role as 
6 did in the exponential distribution, and yvy which is the 
smallest value the random variable can assume. The 
probability density function for the Weibull distribution 


1s: 





The Weibull and Exponential distributions are related by the 
fact that if B=1, v=0, and @=l/a the two are identical. 


Therefore if X has a Weibull distribution with parameters 
(X-v) 
a 





B 
D Qa, and Vv; then r-| | has an Exponential 


] 


distribution with parameter 6=1. Or, X=aY'’+v, where X 
] 


and aY’+v have the same distributions, means, variances, 
r(r-1) 


r 





l 
and higher moments. Now, where E[X’]= , and ala 7 for 


the Exponential distribution, 


1 r 
E|X’\=E a | 2.47a 


J 
al ta, 
_-21r)-iler ja 2.47b 


for the Weibull distribution. Since @ must equal one for 


the similarity to work, paar Z+1}ey, 


Letting v=0, these expressions are used to construct 
the moments used in Equations 2.10 and 2.11. The algebra is 
too involved to show here, but the substitutions are made 
just as they were for the Exponential case in Equations 2.42 


through 2.45. After algebraically reducing the expressions, 


functions for the coefficient of variation (vy), measures of 


skewness (y,) and kurtosis a) are shown below. 
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These expressions show that the measures of skewness and 
kurtosis for a Weibull distribution depend only on the value 
of 8B, the shape parameter of the distribution. 

Once the moment functions were derived, they were 
entered on a spread sheet in Microsoft Excel, where values 
of the coefficient of variation, skewness, and kurtosis 
could be calculated for several different values of the 
distribution parameters. Three plots were constructed: (1) 
skewness versus coefficient of variation, (2) kurtosis 
versus skewness, and (3) kurtosis versus coefficient of 
variation. Figures 2.1, 2.2, and 2.3 show the three plots 
with the curves for each distribution. Notice that the 
Normal and Exponential distributions appear as single points 


on some of the plots. 
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Figure 2.1: Moment Functions 
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Figure 2.2: Moment Functions 
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C. TESTING OF THE MOMENT PLOT METHOD 
1. Random Number Draw 
a. A Graphical Statistical System (AGSS) 

To test whether or not the method would 
discriminate between sample population distributions, random 
numbers drawn from known distributions were analyzed using 
the moment plot method. A Graphical Statistical System 
(AGSS), an APL based software package resident on the Naval 
Postgraduate School's main-frame computer, was used to 
generate the random numbers, and to  calculatewieme 
coefficients of variation, and the measures of skewness and 
kurtosis of the samples analyzed. As part of the analysis 
procedure, AGSS makes histograms, cumulative distribution 
plots, and probability plots for each distribution it )iai@s 
to data. Figures 2.4a, b, and c are examples of these 


plots. 


We'BULL DENSITY FUNCTION, N=68 


0 20 0.25 


i 


RELATIVE FREQUENCY 
0.10 @.15 


005 


g 
B 
8 
8 


Figure 2.4a: Histogram and Weibull Density Function 


22 


WEIBULL CUMULATIVE DISTRIBUTION FUNCTICN, N=68 
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Figure 2.4b: EDF and Weibull CDF 


WHBULL PROBABILITY PLOT, N=68 





Figure 2.4c: Weibull Probability Plot 

Figure 2.4a is a histogram of the sample values 
with the Weibull density function. Comparing the two, shows 
how the frequencies of sample values (represented by the 
bars) differ from the expected frequencies of values from a 
Weibull distribution density function. The empirical 
cumulative distribution function (EDF) is shown in Figure 
2.45 with the Weibull cumulative distribution function 


(CDF). If this sample had a perfect Weibull distribution 
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the empirical plot would lie on top of the Weibull funcuweae 
Figure 2.4c is a Weibull probability plot. It Cconveysmeme 
same information as Figure 2.4b, but is used to more easily 
judge the goodness of fit. If the sample had a perfect 
Weibull distribution, the points plotted would lie on the 


Straight line: To achieve the straight line jie 


probability; y-axis has been plotted as y=In(-In(1- F(x))), and 


the x-axis becomes In(x). Where F (x) is the Weibull 


n 


cumulative distribution function [Reto 
b. Types of Data to be Fitted 

To generate random numbers from a specific 
population distribution, AGSS asks for certain input 
parameters. For Normal, and Lognormal distributions, the 
mean and standard deviation are required. For Weibull 
distributions, the shape and scale parameters are required, 
and for Exponential distributions the mean is required. The 
range and variance of typical flight load data and fatigue 
life data were used to determine the input parameters. Two 
means were chosen to represent the loads data, 500 lbs., and 
5000 lbs. Three means were chosen to represent the fatigue 
life data, 100,000, 1,000,000, and 5,000,000 cycles. 
Standard deviations typical of the data were chosen for 
these means, which were 150, 1500, 30,000, 150,000, and 


750,000, respectively. Five separate random number draws 
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corresponding to these means and standard deviation pairs 
were made from each of the four distributions. 


@. Distribution Fitting 


Once the 20 samples of random numbers were 
generated, each sample was fit to each of the four 
distributions. Figures 2.5a, b, c, d, e, f, g, and h show 
the CDF/EDF plots and probability plots for a Weibull sample 


mrt co the four different distributions. 
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Weibull Sample fit to Lognormal Distribution 
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Figure 2.5e: EDF/CDF Figure 2.5f: Probability 
Weibull Sample fit to Weibull Distribution 
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Figure 2.5g: EDF/CDF Figure 2.5h: Probability 
Weibull Sample fit to Exponential Distribution 


AGSS calculated the sample mean, = standard 
deviation, skewness, and kurtosis for each sample (20 sets 
of values). It also calculated four goodness of fit 
measures for each fit (80 sets of values). 

d. Goodness of Fit Measures 
AGSS calculates the Chi-square, Kolmogorov- 


Smirnov, Cramer-von Mises, and Anderson-Darling goodness of 


fit statistics. The Chi-square goodness of fit test is 
applicable to discrete or continuous, univariate or 
multivariate data, and is the most widely used goodness of 
fit test. There are several variations of the test, all 
based on the following relation, which is a measure of the 
difference between the observed and the expected frequency 
of a value. A graphical representation of this difference 


is Figure 2.4a. 


yay omy) um) Del 


i=] P; 


Here, N is the observed frequency of the ith cell. The 
expected frequency is np, based on the probability density. 
The other three goodness of fit statistics are “EDF 
statistics", based on measuring the vertical difference 
between the empirical distribution function of the sample, 
F(x), and the cumulative distribution function of the 


distribution being tested, F(x). Figure 2.4b shows the 


difference between F (x) and F(x). These vertical 


differences are divided into two classes, the supremum class 
and the quadratic class. 


The supremum statistics are shown in Figure 


2.4b, where D’ is the largest vertical difference when F (x) 


is greater than F(x), and D is the largest vertical 


Ze 


ditterence when F (x) is smaller than F(x). The most well- 


known EDF statistic is the Kolmogorov statistic D. It is 


defined as: D=sup, F.(x)- F(x) 








, and measures the maximum 


difference between F(x), and F(x). {[Ref. 5 and Rei iby 


The quadratic class is given by: 


O= nf {F(x)-F(x)} y(x)aF(x) 2.52 


where y (x) is a weighting function for the s@quaped 
2 

difterence 1 F,(x)- F(x)} . The Cramer-von Mises statistic W’, 

is a special case of the quadratic class where w(x)=1. For 


the Cramer-von Mises statistic, Equation 2.52 specializes 


to: 


w? =n [{E(x)- Fd} d(x) 2.53 


_—o 


The Anderson-Darling statistic A’, is another special case 


of the quadratic statistic where v(x) =[{F(x) H(1- F(x) . FoR 


the Anderson-Darling statistic, Equation 2.53 specializes 


CO: 


,_ t{E,(x)- F(x) f dF (x) 
a J {F(x)}i(1- F(a) lx) 


The quadratic statistics are more powerful than the supremum 
class of EDF statistics, since better use is made of the 
information contained in the whole sample rather than by 
using only the maximum discrepancy [Ref. 6]. Smaller values 
of the four goodness of fit statistics indicate a better fit 
momtne tested distribution. 
2. Results of the Study 
a. Moment Plots 

The following three figures show moments of the 
random samples plotted on the moment function templates. 
There are 20 points in each figure representing the five 


random draws from each of the four distributions. 
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Figure 2.6: Coefficient of Variation vs. Skewness 
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Figure 2.7: Kurtosis vs. Skewness 
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Figure 2.8: Kurtosis vs. Coefficient of Variation 
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The values of skewness and the coefficients of 
variation for the random samples plotted in Figure 2.6 fall 
close to their correct distribution curves. Where the 
Normal and Weibull functions intersect at (.25,0) there are 
three points from the Normal random sample and three points 
from the Weibull random sample. This intersection is the 
only area of ambiguity for this particular plot. Aside from 
this area, the plot indicates the correct distributions of 
the random samples. The random sample values of skewness 
and kurtosis (Figure 2.7) all fall near the curve for a 
Weibull distribution, but the kurtosis is consistently 
underestimated for Lognormal. This plot may not be as 
useful an indicator of the distribution as Figure 2.6. 
Figure 2.8 has three intersection points that could result 
in ambiguous distribution indications. By using all three 
plots, it was hoped that the areas of ambiguity could be 
resolved. 

b. Goodness of Fit Tables 

Goodness of fit statistics were calculated for 
each fit of the random samples and entered in the following 
tables. mon fcance levels are also listed under their 
respective goodness of fit statistic. AGSS specifies a 
Significance level a=0.0l for the goodness of fit 
statistics. It uses as the null hypothesis, H,: “The 


goodness of fit statistic is small enough to indicate a good 


eel 


fit". For each goodness of fit statistic, AGSS calculates 


the corresponding p-value from the sample data. If poae 


the null hypothesis must be rejected (Ref. 7]. Therefore, 
ps90.01 indicates a lack of fit. Tables 2.1 through 2.5 are 
goodness of fit statistics for samples drawn from a Normal 
distribution, then fit to all four distribution Tables 
2.6 through 2.10 are for samples drawn from a Lognormal 
distribution. Tables 2.11 through 2.15 are for samples from 
a Weibull distribution, and Tables 2.16 through 2.20 are for 


samples drawn from an Exponential distribution. 


TABLE 2.1: NORMAL SAMPLE 1 


[Significance ]0.12993_|0 | -+([0.33585 |0. 
[Kolm-Smirn | 0.02433 [0.077617 [0.016414 |0.37174 
[Significance | 0.59467 | 0.000011701 [0.95044 [0 
>0.15 
fAnder-Darl | 0.42392 | infinite [0.34355 [233.6 
0.15. [0 +} >0.15 |<0.01. 















TABLE 2.2: NORMAL SAMPLE 2 


Test \Fit Weibull | Exponential 
Chi-Square 8328 6.0472 | 1611 













5 ; 
0.66595 [0 ~—~*+(0.64194 [oS 
(Significance |0.99351 |0.0000047422 [0.7128 [0 
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TABLE 2.3: NORMAL SAMPLE 3 


Test \Fit 






2.7316 
[Significance |0.95004_|0 +|0.00494 |O. 
0.000092067 [0.95485 [0 
1.8326 
< ; 
1 





Significance 0.01 
0.19491 1.707 0.3155 229.1 
Significance 






TABLE 2.4: NORMAL SAMPLE 4 


Test\Fit Weibull Exponential 
Chi-Square 7.6794 40.694 B00 4144 
Significance |0.36171 |9.267E-7 |0.000069474 |0 — | 






[Significance |0.92973 | 0.025173 [0.041417 [| O 


Kolm-Smirn 0.017169 | 0.046771 0.04403 0.48245 





TABLE 2.5: NORMAL SAMPLE 5 


Test\Fit Weibull Exponential 






(Significance |0.32068 | 0.0013961 |0.0000029004|/0. 
[Significance |0.76541 [0.2145 [0.036249 [0 
fander-Darl _|0.32459 [1.8754 [4.0146 [331.03 






TABLE 2.6: LOGNORMAL SAMPLE 1 


[Significance]0 [0.63735 [0 [0 
(Significance]0—+(|0.74231 [0 [0 
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TABLE 2.7: LOGNORMAL SAMPLE 2 


[Significance|O0_ [0.79192 [0 |o. 




















Significance]O 0 |0.78304 [0 7 


TABLE 2.8: LOGNORMAL SAMPLE 3 


Weibull 
Chi-Square ee reeks 9.989 ISish 6 BZ 1654.4 










[Significance]0 [0.12512 [0 -+|o. 
0.087916 
3,86565-7 [0.34853 |2.0068E-7 [0S 


TABLE 2.9: LOGNORMAL SAMPLE 4 


[Chi-Square [33.105 [4.4023 [121.07 |4021— 

[Significance |0 [0.6224 (0. -+(|o. 

[Significance |0.010641 [0.74173 [0 |o. 
0 















64474 
$0.15 






TABLE 2.10: LOGNORMAL SAMPLE 5 


[Significance]0 [0.71832 [0 +|0. 
[Kolm-smirn | 0.047072 | 0.023323 [0.080214 |0.49845 
[Significance] 0.023791 [0.64815 [0 [0 













Significance | <0.01 
2.5507 | 0.38184 338.83 
Significance 
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TABLE 2.11: WEIBULL SAMPLE 1 


Test \Fit Exponential 
1667.7 
[Significance |0.26697 [0 [0.44025 

0.37167 
47.135 
















Significance <0.01 
0.58359 | 9.7349 0.39178 | 230.33 
Significance <0.01 





TABLE 2.12: WEIBULL 2 


Weibull | Exponential 
Chi-Square 13.324 95.08 9.1728 1490.8 
Significance | 0.1485 4.1633E-17 | 0.32793 
0.02736 | 0.056852 | 0.025149 

0 : 
0 











0.36073 


44251 
44.762 


TABLE 2.13: WEIBULL SAMPLE 3 


Test\Fit Weibull | Exponential 
Chi-Square HD 4. OW 128.99 8.3147 1474.9 
Significance ]0.085657 [0  ~——s—“§“ssd[0.50276 [0 













0.000048364/0.86081 |O 


TABLE 2.14: WEIBULL SAMPLE 4 


Weibull 
Chi-Square 4123.7 3.8461 173.08 66.461 
Significance [0 |0.87073. 10 | 2.4798E-11 












0.017421 
Significance |0 [0.92191 | 6.5556E-7 | 0.0052128 
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TABLE 2.15: WEIBULL SAMPLE 5 


Weibull 
116.44 
[Significance ]0 [0.75813 {0 |1.5068B-10 
[Rolm-Smirn [0.47111 | 0.021057 
[Significance |0 [0.767 | 2.0498E-7 |0.0073019 __ 















TABLE 16: EXPONENTIAL SAMPLE 1 


[Rolm-Smirn [0.16619 | 0.08944 | 0.02144 [0.022324 
5.3839E-25 0.70129 














0.065358 
>0.15 
48.045 0.42406 
0.01 >0.15 








TABLE 2.17: EXPONENTIAL SAMPLE 2 














(Significance|0 +(|0~—~*(| 0.25674 [0.12958 
6.3655 
0.15 





TABLE 2.18: EXPONENTIAL SAMPLE 3 










[Significance ]0 | 1.6945B-14 [0.52516 [0.63775 
[Significance | 1.3516B-20 |0.00006209 [0.92377 |0.92353_ 
44.703 
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TABLE 2.19: EXPONENTIAL SAMPLE 4 


[Significance |0 | 1.3032B-12 [0.43176 |0.34536 
0.070094 | 0.017358 
0.15 
0.01 
















TABLE 2.20: EXPONENTIAL SAMPLE 5 


Test\Fit 










Exponential 
Chi-Square S03 70.603 9.1229 9.5086 
Significance 3.0756E-13 | 0.10427 0.14692 


c. Vote Counting Method 


Ca 
0.16368 0.069234 0.017212 | 0.017594 


The moment plots provide a method of visual 
comparison of the candidate population distribution moments 
with the sample moments, while the goodness of fit 
Statistics provide an analytical comparison that is more 
abstract. The following vote count method combines the two 
procedures to produce a quantitative indication of the 
likely population distribution. 

From the three moment plots (Figures 2.6, 2.7, 
2.8), each of the 60 points provided a "vote" for a 
distribution, based on which distribution moment function 
the point was closest to. If the point fell equidistant 


between two functions, both corresponding distributions 


oy) 


would receive a vote. The goodness of fit statistics 
"voted" for the distribution fit that produced the smallest 
statistic value; i.e., each table above produced four votes. 
The four tables below show the results of this vote counting 
method for each of the known distributions. The first three 
rows are the votes of the moment plots: the next four rows 


are the votes of the goodness of fit statistics. 


TABLE 2.21: RANDOM SAMPLES FROM NORMAL DISTRIBUTION 


TEST \DISTRIBUTION Weibull | Exponential 
S. vs. C. of V. 
Chi-Square =e 
oes ee 
ae 
p26 
















[Kolm.-Smirn. _-‘[ 3 
[Ander.-Dar.—~+«|~—4 


TABLE 2.22: RANDOM SAMPLES FROM LOGNORMAL DISTRIBUTION 


TEST\DISTRIBUTION Weibull | Exponential 
S. vs. C. of V. | 


Kurt. vs. Skew. 


es) 
ae es 
Kurt. vs. C. of v. | 0 | 
ae) 
EA) ae 
Ee ae 
(= OS 
ae 
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Chi-Square 


3 _ _ 
8 | 0 |. 





Vote Totals: 


Pee a Eee 


TABLE 2.23: RANDOM SAMPLES FROM WEIBULL DISTRIBUTION 


Se vs. CC. of V. 
Kurt. vs. Skew. 





TABLE 2.24: RANDOM SAMPLES FROM EXPONENTIAL DISTRIBUTION 
TEST\DISTRIBUTION Weibull | Exponential 
S. vs. C. of V. 


Kurt. vs. Skew. 
Kurt. vs. C. of V. 


ee ee Oe eee ee 
_ Oe ee ee eee 
_ Ua SOS Se ae a ae 

(Chi-Square | 0] a 
i ie (Oe EE ee ee 
Oe Ee | ee 
0 72S eee ee 
ae ee 






Kolm.-Smirn. 


Vote Totals: 





This distribution selection method selected the 
correct distribution in every case but the Exponential. 
This is due to the fact that the moment curves for the 
Weibull distribution intersect the point for the Exponential 
distribution. The Exponential distribution is actually a 
sub-set with specific parameters of the more general Weibull 
distribution. Therefore, Exponential samples’ scattered 
about the Exponential point have a tendency to be closer to 
the Weibull curve than the Exponential point. The selection 
method correctly selected the Normal distribution by a small 
margin over the Weibull distribution. Here, the goodness of 
fit statistics provided more accurate votes than the moment 
plots. 


d. Varying the Sample Size 


The test just described was based on drawing 
Sample sizes of 1000 from a known distribution. The actual 
flight loads and fatigue life data to be analyzed with the 


method have sample sizes around 20. To determine how 


E)y, 


accurate this distribution selection method is with smaller 
sample sizes, the entire test was performed two more times. 
First, random samples of 100 were drawn from known 
distributions, analyzed with AGSS, the moments were plotted, 
goodness of fit statistics calculated, and the votes 
tallied. Then random samples of 20 were drawn, and the 
process was repeated. Figures 2.9a, b, c, and d show 
EDF/CDF plots and probability plots for Weibull samples fit 
to a Weibull distribution. Comparing these four Figures to 
Figures 2.5e and 2.5f gives a preliminary indication of how 


sample size affects the goodness of fit. 


WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=100 WEIBULL PROBABILITY PLOT, N=100 


10 
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0.8 
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0.2 
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800 BOC 1000 100 


Figure 2.9a: EDF/CDF Figure 2.9b: Probability 
Random Weibull Sample of 100 fit to Weibull Distribution 


40 


WEIBULL CUMLLATIVE DISTRIBUTION FUNCTION, N=20 WEIBULL PROBASILTY PLOT, Nw=20 


08 
~g 
tf 


0.6 
PERCENT) 
my 


0.4 


CUMULATIVE PROBABILITY 


0.2 





200 400 6co aco 400 
SAMPLES FROM WEIBULL OfSTRIBUTION SAMPLES FROM WEIBULL DISTRIBUTION 


Figure 2.9c: EDF/CDF Figure 2.9d: Probability 
Random Weibull Sample of 20 fit to Weibull Distribution 


Only the total vote count tables for each distribution are 
shown below. Tables 2.25 through 2.28 are for sample sizes 


of 100, Tables 2.29 trough 2.32 are for sample sizes of 20. 


TABLE 2.25: SAMPLE SIZE OF 100 FROM NORMAL DISTRIBUTION 


Moment Plots [| 11 | 2 | 7 | 0 
(Goodness of fit_|. 7. |3. | ~6 ~| 0 
vote Totals: [| 18 | 5 | 13. 0 






TABLE 2.26: SAMPLE SIZE OF 100 FROM SUNS DISTRIBUTION 


Meericts [| 5 |. 4 | 6 | o 
(Goodness of fit | 0 | 16 | 4 ~~| 0 
[Vote Totals: | 5 | 20.~+| 10 ~+4~ SOS 






TABLE 2.27: SAMPLE SIZE OF 100 FROM WEIBULL DISTRIBUTION 


[Moment Plots | 4 | 1 | 10. 0 
[Goodness of fit | 9 | 0 | a1 [| 0 
[Vote Totals: | 13. [.1~~+|~21+~*+| 0S 
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TABLE 2.28: SAMPLE SIZE OF 100 FROM EXPONENTIAL 












DISTRIBUTION 
TEST\DISTRIBUTION | Normal | Lognormal | Weibull | Exponential 


Moment Plots a 
Goodness of fit [|.o0 | 0 | 16 | 4 
[vote Totals: | 0 | 0 | 30 | 41 









TABLE 2.29: SAMPLE SIZE OF 20 FROM NORMAL DISTRIBUTION 


[Moment Plots | 5 | 2 | 10 | 0 
[Goodness of fit | 5 | 11.,-=~) ~4~Y| ~O Cs 
[vote Totals: | 10. | 13. | 14 ~~) 0 1 











TABLE 2.30: SAMPLE SIZE OF 20 FROM LOGNORMAL DISTRIBUTION 


[Moment Plots | 3 | 6 | 6 | 0 
[Goodness of fit [| 2. | 11 ~=)~+4~+|~ 0 
[vote Totals: | 5 | 17.~+«Y~-a—~<Y(C™S~<“—~S SCS 











TABLE 2.31: SAMPLE SIZE OF 20 FROM WEIBULL DISTRIBUTION 


(Moment Plots [5 | 3. ~+| 13] 0 
[Goodness of fit | 5 | 3. ~+| 9 ~+|~ ~~ 0~S=S—~s™S 
[vote Totals: [10] 6 | 22 | 0 






TABLE 2.32: SAMPLE SIZE OF 20 FROM EXPONENTIAL DISTRIBUTION 


[Moment Plots | 0 | 0 | 15 | 2 
[Goodness of fit | 0. | 3 ~|~7~+)~ SOC‘ 
[vote Totals: | 0 | 3 | 22 | 7 | 







The vote count method was just as effective with 
sample sizes of 100 as with 1000. Lognormal and Weibull 
distributions were correctly selected by a wide margin. 


Again, the Normal distribution was correctly selected by a 
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small margin over the Weibull distribution with the moment 
plots indicating with greater accuracy, while Lognormal was 
also correctly selected, the goodness of fit tests were more 
accurate. Weibull was correctly selected with the moment 
plots indicating most accurately, but both methods missed on 
the Exponential. With sample sizes of 20 the method loses 
considerable accuracy. For the known Normal samples, moment 
plots indicated Weibull, goodness of fit statistics 
indicated Lognormal. However, the total votes for Normal 
was close to the totals for Lognormal and Weibull. 
Lognormal was correctly selected by a 3:2 margin, with the 
goodness of fit statistics providing the most accurate 
votes. Weibull was correctly selected by both components of 
the method, with a 2:1 margin. Based on the three test 
results of varying sample sizes, sometimes the moment plot 
is better, sometimes the goodness of fit tests are better. 
Combining the two, the distribution selection method is able 
to correctly determine whether a sample is from a Weibull or 
Lognormal distribution in most cases. When a sample has a 
coefficient of variation between 0.2 and 0.5, the moment 
plots method cannot differentiate between Normal and 
Weibull. The reason for this is that intersections occur 
for the Weibull and Normal distribution moment functions 
within this range (see Figure 2.4 and Figure 2.6). It would 


be reasonable to assume that if two distributions receive 
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Similar vote counts for a sample being analyzed, either 
distribution could be used to model the sample's population. 
Another alternative in case of a tie would be to tally the 
votes again, allowing each moment point and goodness of fit 
statistic to vote for one of the tied distributions only, 


after eliminating the others. 
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III. APPLICATIONS OF THE MOMENT METHOD 

A. HELICOPTER FLIGHT LOADS DATA 

The first application of the moment/goodness of fit vote 
count method was to flight loads measured on the main rotor 
forward longitudinal stationary star (MRFLSS) of the SH-60 
Sikorsky helicopter. This is a fatigue critical component. 
Data was provided from two separate substantiation flights 
flown by United Technologies Corporation, Sikorsky Aircraft 
Division. Load meaSurements were recorded for a variety of 
maneuvers at a variety of airspeeds, collective settings, 
and "g" loadings. The maneuver that caused the most fatigue 
damaging loads was the symmetric pull-out. Table 3.1 lists 
the two flights at two different gross weights, and the runs 
Mage Curing each flight. Essential parameters of each run 
are also listed. Parameter abbreviations are: SP is 
symmetric pullout, 155K or 124K is the airspeed in knots, FC 
is fixed collective, TC is top collective, -25% is the 
collective set at 75% of TC, and G's are the accelerations 
the maneuver produced in terms of the acceleration of 


gravity. 
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TABLE 3.1: FLIGHT/RUN PARAMETERS 
Flight 37: Flight 43: 
16,500 gross wt. 20,800 gross wt. 
_ MANEUVER 


_ RUN | MANEUVER |_CRUN 
14 SYM PULL (SP), 124K SP, 124K, 100% 
FIX COLL (FC), a2mcG TORQUE (100%), Ieee 

15 (SP, 124k perecrcrcc SP, 124K, 100%, 1.75G 
SP, 124K, FC, 2.8G SP, 124K, 100%, 2.0G 
SP, 247, Ss 0c SP, 124K, 100%, 2age 
SP, 124K -25%, 2.9G SP, 124K, 100%, 2.1G 
SP, 124K -25%, 2.8G SP, 124K, 1008)020e 
| 22 |SP, 124K, -25%, 3.0G SP, 124K, 100%, 1.75G 
SP, 124K, -25%, 2.96 SP, 124K, 100%, 2.1G 


[24  |SP, 124K, -25%, 2.9G SP, 124K, 100%, 2m 


25 | SP, 124K, TOP COLL 81 | SP, 155K, 100%, 1.5G 
(TG). 26G 


| 28 ESP eS 5K = 258m oG SP, 155K, 1003) 
29 SP, 155K, =258) ose SP, 155K, 100%) 92a 
| 30 ESP, ISSK — 258 ore SP, 155K, 105% /) ee 
| 31 SSP AISSKy a2 oe) oe SP, 155K, 105%, 2am 
| | ———<“C«~*™*C‘C:*CS*YCSC* SP, TSK, LOC ey eee 
po 8 SP, 155K, 100%, 1.856 
po 8B SP, 155K, 100%, 2.36 














1. Loads Data Processing 
From a statistical standpoint, it would be ideal to 
have several runs made for each set of parameters; rather 
than a different set of parameters for each run. HoOwewem® 
this is not deemed economically feasible with substantiation 
flights at the present time; therefore, the flight loads 
data had to be pooled prior to any statistical analysis. 
a. Sorting and Pooling Data 
Load versus time plots were made of each run. 
Flight 43, which had the higher gross weight, produced 


higher loads, even though the "g" loading was slightly less 
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[Eman flight 37. Within the two flights, the runs at 155 
kts. produced slightly higher loads than those at 124 kts. 


Figures 3.1 and 3.2 are examples of the load versus time 


rots. Figure 3.1 is at 124 kts.: Figure 3.2 is at 155 
kts. 
4500 
4000 ! 
3500 ‘| 
q 3000 \ | 
= 2500 | 
” 
8 2000 
5 1500 \ 
1000 i 
geo AI pn Vay a any h 
NA Pm eee 
0 , 
SNreenmenesgsezsergasrsseanesyg & B 
eee eo NS KH NRK OK UO 
vo) S = “ © wy SS ©. - - =- = Saye Sey So a Se 
© 
Time (sec.) 
Figure 3.1: 124 kts. Load vs. Time 
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Figure 3.2: 155 kts. Load vs. Time 


Initially four pools were established, one for each fivem: 
speed within the two flights. 
b. Selecting Peak Loads 

Fatigue life calculations only consider the 
maximum and minimum loads in a given load cycle. Because 
these loads were sampled every 0.17 seconds, many of the 
loads recorded were intermediate loads en route to a maximum 
or minimum. The loads of each run were processed through a 
computer algorithm to pick out only the maximum lodadS?agem 
"peaks", which are the loads of principal interest. Tee 
needed, the same analysis could be carried out to find the 
Minimum loads in each excursion, or what are called the 


sve ey S 
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c. Selecting Cut-off Loads 

In addition to helping pool the data, the load 
versus time plots indicated that the recording instruments 
were turned on just before the maneuver was initiated, and 
turned off just after the maneuver was completed. Loads 
encountered in straight and level flight are not fatigue 
damaging, and are not necessarily from the same population 
as the fatigue damaging loads produced during the maneuver. 
A minimum cut-off load was determined for each pool of data. 
Loads less than the cut-off load were not considered in the 
analysis. The cut-off load was determined by plotting the 
empirical cumulative distribution function for the loads of 
each run. Each of these plots produced a significant cusp, 
which indicated the possibility of two populations. One for 
the straight and level flight loads and one for the maneuver 
loads. The peak of the cusp is taken as the cut-off load. 
(Ref. 8] Figures 3.3 and 3.4 are typical empirical 
Cumulative distribution function plots. Figure 3.3 is a 
plot of all the recorded loads. Figure 3.4 is a plot of 


only the maximum loads. 
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Figure 3.3: ECDF for all Loads of Flight 43 Run 88. 
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Figure 3.4: ECDF for Peak Loads of Flight 43 Run 88 
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The ECDF's for the flight 43 runs indicated a cut-off load 
of 750 lbs. for both flight speeds. The ECDF's for the 
flight 37 runs indicated cut-off loads of 500 lbs. for runs 
flown at 124 kts. and 1400 lbs. for runs flown at 155 kts. 
2. Selecting the Best Distribution 

Once the peak loads above the cut-off were selected 
from each run, all the runs were analyzed with AGSS to 
obtain the moment values and goodness of fit statistics. 

a. Moment Plots 

Three moment plots were made for each of the two 

Ieergncs. Figures 3.5, 3.6, and 3.7 are for Flight 37, and 
fmeces 3.8, 3.9, and 3.10 are for Flight 43. Each data 


point on the plots represents a run from that particular 


flight. 
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Figure 3.5: Flight 37 Coefficient of Variation vs. Skewness 
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Figure 3.9: Flight 43 Kurtosis vs. Skewness 
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Figure 3.10: Flight 43 Kurtosis vs. Coefficient of 
Variation 


b. Goodness of Fit Results 


Tables 3.2 through 3.15 contain the goodness of 
fit statistics for flight 37. Tables 3.16 through 3.32 are 
for flight 43. Just as each run was represented by a point 
on the moment plots, each goodness of fit table represents 


the indicated run. 
Table 3.2: Flight 37, Run 14 
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TABLE 3.3: FLIGHT 37, RUN 15 


[Significance |0.76353 | 0.85947 
fAnder-Darl [0.6294 | 0.71922 
50.15 














TABLE 3.4: FLIGHT 37, RUN 16 


4.1774 
0.040967 
0.15847 












0.15 
0.65443 
>0.15 





TABLE 3.5: FLIGHT 37, RUN 17 


Test \ Fit Weibull | Exponential 














[Cramer-V M__|0.14312 [0.16876 [0.1522 [0.19064 _| 


Ander-Darl 0.91033 i002 7 0.94135 ie 2io7 
Significance 


TABLE 3.6: FLIGHT 37, RUN 19 


0.15 
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TABLE 3.7: FLIGHT 37, RUN 21 


Normal 
[chi-Square [| | ~~ —*+dY~SCS TS 
[Significance |_—~+|——~«dY S739” 










TABLE 3.8: FLIGHT 37, RUN 22 


[Chi-square | ~*~. SSSSC~idSC—CSCS~C~S~S 233 
[Significance] +(|—~—~—~S~SYS—CSCS~CS~ STC” 
Kolm-Smirn 
















TABLE 3.9: FLIGHT 37, RUN 23 


[chi-square | +. ——SC<~SSC<CS~S TTC 
[Significance | -+(|~—~——~—S«(Y Ss 00052727 















TABLE 3.10: FLIGHT 37, RUN 24 


7 
(Chi-square ]2.9906 | ——~—~+#2.1593 2.0009 
[Significance |0.08375 |__| 0.14171 0.1572 
[Kolm-Smirn | 0.25262 
[Significance | 0.29301 [0.74016 
[Cramer-V M__] 0.19108 












Test \ Fit 





B15, 


TABLE 3.11: FLIGHT 37, RUN 25 


12.977 8.783 
0.0015209 0.066747 
0.15802 0.20625 
0.15 0.15 

0.025 







TABLE 3.12: FLIGHT 37, RUN 28 


1.6221 [2.6723 [2.1401 [65.469 
0.20279 



















0.034765 
>0.15 
>0.15 >0.15 


TABLE 3.13: FLIGHT 37, RUN 29 


Kolm-Smirn | 0.11639 | 0.15903 


















Significance 


TABLE 3.14: FLIGHT 37, RUN 30 


Test_\ Fit 
Chi-Square 1FS03 2 1.1959 51.488 











Ander-Darl 0.50851 |0.88421 0235795 6.7205 
Significance 







oie | 


TABLE 3.15: FLIGHT 37, RUN 31 


0.22167 0.41391 

0.1713 
: 0 
> 





15717 





TABLE 3.16: FLIGHT 43, RUN 69 






Weibull | Exponential 
|Chi-Square | | 
Significance 






ee 








TABLE 3.17: FLIGHT 43, RUN 70 


[chi-square [| SO COCSCSCS~*SdSCSS 
[Significance | —~«|~—SSCS< ] SSSSCidS SS 











TABLE 3.18: FLIGHT 43, RUN 71 


[chi-square |) SOS TTCOCrYS—CSCSC 
[Significance | ——S=«Y.SSCSCirSCSCS~S~idSSC(‘<CSCS~S~S~S~S 
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TABLE 3.19: FLIGHT 43, RUN 72 


Test_\ Fit [Normal | Lognormal |Weibull 
Chi-square | | *+4|+d|SSSSCSCSCS 
Significance | -+| ~~ ~*+|—S«<Y‘~CSCSCSCS<«7<;«< 















0.73248 |0.6453 0.67992 | 0.47017 
0.099269 |0.11644 0.10509 | 0.25719 
On 15 AO ILS >0.15 SO) Sis 


Significance 


TABLE 3.20: FLIGHT 43, RUN 73 


[chi-square | | SY SSSCCCSdS—<~;~CSCSCSCSCSS 
[Significance | ~*<«| ~~. ~SCS~S/ SCStCSCSY 

















TABLE 3.21: FLIGHT 43, RUN 74 


(chi-square | | +4. SSS 
(Significance | ——~«| ~~ SSSC~dY'SCOSCdCOTTCSOSCS 
50.15 50.15 





TABLE 3.22: FLIGHT 43, RUN 75 


[chi-square | | || | 
[Significance | | | 
yd gs SO. 15 yeris yoru 













51, 


TABLE 3.23: FLIGHT 43, RUN 76 


[chi-square {S| t—“‘i ee 
[Significance | +(| ~~. ~SC<CS~<~dX SSSC<C~;7;7<CS*# 


















TABLE 3.24: FLIGHT 43, RUN 77 


Lognormal 
[chi-square | | SSS—CSs—<‘dS—SC 
[Significance | *+| ~~ ~*<Y|.SdY—aSCSCSCS<;«;7«<C 
[Significance | 0.85407 | 0.41875 







Significance 
TAnder-Darl |0.41626 |0.58914 0.48493 [0.9025 
Significance |>0.15 




















TABLE 3.25: FLIGHT 43, RUN 81 


Significance | >0.15 





TABLE 3.26: FLIGHT 43, RUN 82 


Weibull 
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TABLE 3.27: FLIGHT 43, RUN 83 


2.7166 
0.25709 
0.20433 0.16294 





Significance 0.15 
Ander-Darl 0.88484 0.49288 0.67285 3.1335 
Significance 0.025 









TABLE 3.28: FLIGHT 43, RUN 84 


[chi-square [| |. | +. SSCS 
[Significance|—~*+(| ~~ SY|SSCSCS~SYStC‘<CS;<;<S;<7;<;<;<;<C 


















TABLE 3.29: FLIGHT 43, RUN 85 


Weibull | Exponential 
chi-square [| | SSYSSSTCSCSYS—“CSCSTSCST 
Significance | | 


Saar 
[Significance |>0.15 0 [>0.15 [90-15 [0-15 
Significance 









Test \ Fit 











TABLE 3.30: FLIGHT 43, RUN 86 


Test \ Fit Weibull | Exponential 
Chi-Square | ——~+Y~S<«C<SC«‘<S 6.0039 








ia Ee 
Significance | ee ez 
0.20364 | 0.27296 .1564 0.5081 







0 
[Significance |0.60712 [0.24786 | 0.86324 
0.098102 0.049068 
50.15 
Significance [>0.15 
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TABLE 3.31: FLIGHT 43, RUN 87 


Weibull [Exponential 



















0.20776 
0.15245 


TABLE 3.32: FLIGHT 43, RUN 88 


[chi-Square | | SSO] 
[Significance] —~«| SSCS 






Co." Vote Takiy, 
Tables 3.33 and 3.34 contain the vote counts for 


flights 37 and 43 respectively. 


TABLE 3.33: FLIGHT 37 VOTE COUNTS 


[Moment Plots | 13 | 0 | 30 | 0 
[Goodness of fit_| 18 | 16 | 12. | 6 
[vote Totals: | 31 [| 16 | 42 | 6 






TABLE 3.34: FLIGHT 43 VOTE COUNTS 


TEST\DISTRIBUTION Weibull | Exponential 






Moment Plots [iis | 0 | 30 | .— 
Goodness of fit 
Vote Totals: 
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The goodness of fit statistics chose a Normal 
distribution for both flights by a small margin, while the 
moment plotting method chose Weibull for both flights by a 
large margin. The total votes indicate that the flight 37 
load samples are from a Weibull distribution, while the load 
samples from flight 47 are from a Normal distribution. The 
methods having eliminated all but Normal and Weibull, the 
votes were tallied again, allowing the points and goodness 
of fit statistics to vote solely for Normal or Weibull 
distributions. The following tables show the resulting vote 
totals. 


TABLE 3.35: FLIGHT 37 VOTE COUNTS 


TABLE 3.36: FLIGHT 43 VOTE COUNTS 










These tables strengthen the case for a Weibull distribution 
for flight 37, and raise doubts about flight 43 having a 
Normal Distribution. It seems unlikely that simply changing 
the gross weight of the aircraft would change the load 
Sestribution. If curve fits were made of the data, it is 


Clear that plots on Figure 3.5 would indicate Weibull-like 
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behavior with a positive slope. On Figure 3.6 a CUrvVemeae 
indicates Weibull-like rather than Lognormal; thus the 
moment method strongly indicates Weibull for flight 37. 
Similar arguments point to Weibull for flight 43 also; thus 
a curve fit approach would confirm the moment methods vote 
results. 
C. A-6 COUNTING ACCELEROMETER DATA 
1. Explanation of Data 

The third application of the moment plot/goodness of 
fit vote counting method was to counting accelerometer 
measurements from a pool of 103, A-6 aircraft. The 
accelerometers count the number of exceedences of 4, 5, 6, 
and 7 "g" accelerations. They are used to track the fatigue 
lives of each aircraft by recording the fatigue damaging 
loads each aircraft experiences. The recorded loads are 
then converted, through methods of fatigue analysis, into a 
percentage of fatigue life expended (FLE). Once an aircraft 
reaches an FLE of 0.67, it has consumed 67% of its fatigue 
life and is restricted to lower flight accelerations. [Ref. 
9] The pool of 103 were extracted from a population of 351 
A-6 aircraft in the Navy's inventory that possessed usable 
data. The 103 selected were all the flight unrestricted 
aircraft in the population. The values of the random 
variable were the number of exceedences of a given g load 


per aircraft per 1000 flight hours. The data were converted 
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from exceedence data to occurrences, where 7 g's do not also 
count as occurrences of 4, 5, and 6 g's. 
2. Moment Plots 

AGSS was used to compute sample values of the 
coefficient of variation, skewness, and kurtosis for each g 
loading. This application of the moment plot method used a 
different region of the template. The number of g 
occurrences varied greatly between aircraft, particularly at 
the 7 g level, where only a few aircraft had any occurrences 
ae all. This produced large standard deviations. Sample 
coefficients of variation ranged from 0.6 to 2.47, compared 
to coefficients of variation for the helicopter flight loads 
which ranged from 0.1 to 0.9. 

Figures 3.11 through 3.13 are the moment plots. 
Each plot has four points representing the four different g 
Joads. Each point is calculated from the data of all 103 


aircraft. 
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Figure 3.11: Skewness vs. Coefficient of Variation 
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Figure 3.12: Kurtosis vs. Skewness 
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Figure 3.13: Kurtosis vs. Coefficient of Variation 


Lie 


Goodness of Fit Results 


The goodness of fit statistics AGSS calculated are 


Shown in tables 3.37 through 3.40. Each g loading has a set 


of 


values. This data was only fit to the three 


Ggseributions indicated. 


TABLE 3.37: 4G OCCURRENCES 


eibull 
"5232 










W 
5 
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TABLE 3.38: 5G OCCURRENCES 
















TABLE 3.39: 6G OCCURRENCES 
















TABLE 3.40: 7G OCCURRENCES 


[Chi-Square [39.264 | SCSS SSS 
[Significance |3.70188-10| sd. SCS 












3.6886 
0.01 






4. Vote Tally 
Tables 3.41 through 3.44 contain the vote count 


results from the accelerometer data. 
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TABLE 3.41: 4G VOTES COUNTS 


[Weibull 
Moment Plots 
oo 


Zn 
Goodness of fit | 0 | 4 
[vote Totals: | 0 | 6 










TABLE 3.42: 5G VOTE COUNTS 





TABLE 3.43: 6G VOTE COUNTS 


Moment Plots 
Goodness of fit 
Vote Totals: 





TABLE 3.44: 7G VOTE COUNTS 


Moment Plots 


Goodness of fit 
Vote Totals: 





The distribution selection method chose a Lognormal 
distribution for the 4g and 5 g occurrences, and a Weibull 
Meetribution for the 6 g and 7 g occurrences. From the 
earlier random number tests of the method, it was determined 
that the method was most accurate at selecting Lognormal and 
Weibull distributions. That does not mean that the method 
is biased toward those distributions, but rather, when the 


random sample was from a Lognormal distribution, it chose 
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that distribution by a wide margin. The same was true for a 
random sample from a Weibull distribution. 
D. AGARD A-7 COUNTING ACCELEROMETER DATA 
1. Explanation of data 
The following analysis was performed on data taken 
from AGARD Conference Proceedings 506 [Ref. 10]. It was 


chosen because it contained more data for the high g loads, 


than the A-6 application above. Data was available on 40 
aircraft. It was in terms of exceedences and had been 
standardized to exceedences per 1000 hours. The moment 


plot/goodness of fit statistics vote count method was 
performed on the data twice. Once on the exceedence data, 
and again on the data after it had been converted to 
occurrences. The object was to see if the distributions 
were affected by the form of the counts, 1.e., exceedences 
or occurrences. 
2. Moment Plots 

Figures 3.14 through 3.16 are the moment plots for 

the data in exceedence form. Figures 3.17 through 3.19 are 


the moment plots for the data in occurrence form. 
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Figure 3.14: Skew. vs. Coef. of Var. for Exceedence Data 
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Figure 3.15: Kurt. vs. Skew. for Exceedence Data 
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Figure 3.16: Kurt. vs. Coef. of Var. for Exceedence Data 
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Figure 3.17: Skew. vs. Coef. of Var. for Occurrence Data 
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Figure 3.19: 
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for Occurrence Data 


~——#®—— Normal 
——tG-—— Lognormal 
——*—— Weibull 


Exponential 


5g Cor. 


6g Cor. 
7g Cor. 


8g Cor. 





1.5 


for Occurrence Data 


These six figures show that the moment plots are very 
Similar between the exceedence and occurrence data. Curve 
fits of the data in Figures 3.14 and 3.17 would indicate non 
Normal and non Exponential. 
3. Goodness of Fit Results 

Goodness of fit statistics were calculated for the 
data both in exceedence form and in occurrence form. Tables 
3.45 through 3.48 are for the exceedence data. Tables 3.49 


through 3.52 are for the occurrence data. 


TABLE 3.45: 5G EXCEEDENCES 


Test \ Fit Exponential 
Chi-Square 1 HOae C2 eo TORS 62 54.425 
Significance | 0.0053284 | 1.4318E-5 | 5.035E-3 | 9.1051E-12 


















TABLE 3.46: 6G EXCEEDENCES 


Exponential 
54.609 

8. 3185-12 
0.32351 


0.00037485 
1.4046 
0.01 
6.8894 
1 0.01 
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TABLE 3.47: 7G EXCEEDENCES 


[chi-Square_[ ‘| 11.632, ~+| ~~~ ~—~*<~' 3-914 
[Significance | ___-+|5.62E-4_ | ~~~~*«:'[0.047884 
1. 158E-4 0.014466 

< 















1.2139 0.74737 | 0.73626 
Significance <0.025 


Ge 0Zo> 


AI|~I}ATER 
Ole | O 
e FWO le TWN 
OO} | OO} W 
Holla 
\O 


TABLE 3.48: 8G EXCEEDENCES 





TABLE 3.49: 5G OCCURRENCES 


Test \ Fit Weibull | Exponential 
Chi-Square 36 7408) esac aye: 252255 43.361 
Significance | 0.074158 5.1235E-4 | 0.32934 3.8388E-10 














Significance 0.15 





TABLE 3.50: 6G OCCURRENCES 


11.499 53.55 
0.20694 
(Significance | 0.059685 185-4 

0 : 
< 


8 
.39932  |0.33598 | 0.36305 1.3617 

< 

6 






Significance 0.01 
fAnder-Darl _|2.5144 [1.798 [2.0962 
0.05 





i 


TABLE 3.51: 7G OCCURRENCES 


[chi-Square_| [36.482 [5.3093 [3.0649 
[Significance| | 1.5407E-9 [0.021213 [0.080001 










<0. 01 
Ander-Darl 


TABLE 3.52: 8G OCCURRENCES 


Test \ Fit Weibull Exponential 


















[chi-Square [| —~—S~dY—S~<CS;73 P*?)S*S*SS*«~*; CCC 
[Significance |——~«~|~S~<CS~SStC<CS~S«~d 


4. Vote Tally 
Tables 3.53, 54, 55, and 56 show the distributions 
selected for the exceedence data. Tables 3.57, 58, 59, and 
60 show the distributions selected for the occurrence data. 
Both have eliminated Normal and Exponential distributions 


based on curve fits of the data in Figures 3.14 and 3.17 


TABLE 3.53: VOTES FOR 5G EXCEEDENCE 











TEST\DISTRIBUTION 
[Moment Plots | 0 | 1 | 2 | 0 
[Goodness of fit [| 0 | 0 | 4 | 0 
[vote Totals: | 0 | 1 | 6 | 0 
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TABLE 3.54: VOTES FOR 6G EXCEEDENCE 


TEST\DISTRIBUTION | Normal 
Moment Plots 


'Goodness of fit 
Vote Totals: 





TABLE 3.55: VOTES FOR 7G EXCEEDENCE 


TEST\DISTRIBUTION Weibull | Exponential 






Moment Plots ee a ee eee 
Goodness of fit [| 0 | 0 | 3. | 0 
Vote Totals: et a a a 





TABLE 3.56: VOTES FOR 8G EXCEEDENCE 


[Moment Plots | 0 | 1 | 2 | 0 
[Goodness of fit | 0 | 0 | 3 | 0 
[vote Totals: |0 | 1 ~| 5 | 0 






TABLE 3.57: VOTES FOR 5G OCCURRENCES 


Moment Plots | 0 | 1 | 2 | 0 
[Goodness of fit | 0 | 0 | 4 | 0 
[vote Totals: |0 |. 1. | 6 | 0 






TABLE 3.58: VOTES FOR 6G OCCURRENCES 


Moment Plots | 0 | 2 | 1 | 0 
[Goodness of fit | 0 | 4 | 0 | 0 
[vote Totals: [| 0 ]|.6 | 1 | 0 






TABLE 3.59: VOTES FOR 7G OCCURRENCES 


Normal 
Moment Plots.| 0 | 1 | 2 | 0 | 
[Goodness of fit [0 | 0 | 3 | 0 
[vote Totals: | o | 1 | 5 | 0 










ae 


TABLE 3.60: VOTES FOR 8G OCCURRENCES 
TEST\DISTRIBUTION 


Moment Plots 


Goodness of fit 


Vote Totals: 





The 5g, 79g, and 8g data are consistently Weibull 
regardless of whether the data is in exceedence or 
occurrence form. The 6g data are Lognormal in both cases, 
which is intuitively disconcerting. It was expected that 
all the loads would be described by the same distribution 
type. There was some consistency in behavior among g levels 
in A-6 and A-7 data. The A-6 data was Lognormal for the two 
lower g levels, 4g and 5g; and Weibull for the two higher g 
levels, 6g and 7g. The A-7 data was also Weibull for the 
two higher g levels, 7g and 8g; but the two lower g levels, 
5g and 6g, were split between Weibull and Lognormal 
respectively. 

B. SPECIMEN FATIGUE LIFE TESTS 

The second application of the moment plot/goodness of 
fit vote count method was on data generated by fatigue life 
testing of Aluminum 7075-T6 samples. 

1. Explanation of Data 

The fatigue testing was originally performed to 
investigate the affects of mean strain on fatigue lives. 
Twenty samples were first tested at zero mean strain with an 


oscillating strain amplitude of 0.007 in/in. The fatigue 
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lives of each sample were recorded in cycles to failure. 
The strain amplitude was reduced to 0.005 in/in and 20 more 
Samples were tested and their cycles to failure recorded. 
This process was repeated two more times for strain 
amplitudes of 0.003 in/in, and 0.0025 in/in. At this point 
80 lives had been recorded in terms of cycles to failure. 
Next, the mean strain was increased to 0.03 in/in and the 
testing at the four strain amplitudes described above was 
repeated producing 80 more fatigue lives. Two more mean 
strains were tested in this manner: 0.063 in/in and 0.100 
infin. A total of 16 combinations of mean strain and strain 
amplitude were tested, with 20 samples for each combination 
producing 320 fatigue lives. [Ref. 11] 
2. Moment Plots 

The AGSS analysis was made on each of the 16 sets of 
lives, producing values for the coefficient of variation, 
skewness, kurtosis, and goodness of fit statistics. All 16 
tests are shown together on Figures 3.20, 3.21, and 3.22. 
Each point represents 20 samples tested at a _ given 
combination of mean strain and strain amplitude. MThe trend 
in Figure 3.20 is too close to differentiate Weibull from 
Normal, but the negative values of skewness in Figure 3.21 
and the trend in Figure 3.22 eliminate Lognormal as a 


candidate. 
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Figure 3.20: Specimen Skewness vs. Coefficient of Variation 
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Figure 3.21: Specimen Kurtosis vs. Skewness 
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Figure 3.22: Specimen Kurtosis vs. Coefficient of Variation 


Moment plots were also made for two different sub 
poolings of this data to see if mean strain or strain 
amplitude had an affect on the distribution of lives. Four 
pools, each consisting of tests at the same mean strain and 
varying strain amplitudes were plotted first. Then four 
pools, each consisting of tests at the same strain amplitude 
and varying mean strains were plotted. No noticeable trends 
Or Groupings resulted from sub pooling the data; therefore, 
the plots resulting from a single pool of all the data were 
used in the vote counting. 

3. Goodness of Fit Results 
Tables 3.60, through 3.75 contain the goodness of 


fit statistics. Each table represents a test of 20 samples. 
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Table names correspond to variable names assigned in the 


fatigue tests. 


TABLE 3.60: MEAN STRAIN = 0.0, STRAIN AMPLITUDE = .007 


[Chi-Square__[3.5466. | ~——~*+(| 3.5887 [20.827 
(Significance | 0.059668 | ~~—+| 0.058173 | 0.000030036 
















TABLE 3.61: MEAN STRAIN = 0.0, STRAIN AMPLITUDE = .005 
















TABLE 3.62: MEAN STRAIN = 0.0, STRAIN AMPLITUDE = .003 
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TABLE 3.63: MEAN STRAIN = 0.0, STRAIN AMPLITUDE = .0025 


[chi-Square [7.2004 | ~—+|6.93. 7) SS 
[Significance |0.0072559|_~~—~+([0.0084761| SOS 
18653 
48968 










OO}; O]|O;O 





TABLE 3.64: MEAN STRAIN = 0.03, STRAIN AMPLITUDE = .007 


0.46384 
0.49584 
Kolm-Smirn 0.090969 |0.074902 
0 0 
; 0 0 














99643 99987 


TABLE 3.65: MEAN STRAIN = 0.03, STRAIN AMPLITUDE = .005 


Test _\ Fit Exponential 
Chi-Square 0.73182 1.6243 0.6278 NOR 22 







[Significance }>0.15 _[>0.15 [>0.15 [<0.01 | 













TABLE 3.66: MEAN STRAIN = 0.03, STRAIN AMPLITUDE = .003 


Test \ Fit Weibull | Exponential 





[Chi-Square [1.6331 [1.568 [1.3487 |S 
[Significance |0.44194 |0.45656 [0.50948 |S 
Kolm-Smirn 
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TABLE 3.67: MEAN STRAIN = 0.03, STRAIN AMPLITUDE = .0025 


Test \ Fit [Normal |Lognormal | Weibull 
Chi-Square 202 oS pe es Ey, 2.3697 1.1322 








0.15957 







TABLE 3.68: MEAN STRAIN = 0.063, STRAIN AMPLITUDE = .007 











TABLE 3.69: MEAN STRAIN = 0.063, STRAIN AMPLITUDE = .005 











TABLE 3.70: MEAN STRAIN = 0.063, STRAIN AMPLITUDE = .003 
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Test \ Fit 
0.79929 
0.074685 
0.15 70.15 














TABLE 3.71: MEAN STRAIN = 0.063, STRAIN AMPLITUDE = .0025 






TABLE 3.72: MEAN STRAIN = 0.100, STRAIN AMPLITUDE = .007 


[Significance | 0.57106 | 0.78355 
Kolm-Smirn 






>0.15 0.15 






TABLE 3.73: MEAN STRAIN = 0.100, STRAIN AMPLITUDE = .005 


0.056647 
0.81188 
0.10215 











0.15 


TABLE 3.74: MEAN STRAIN = 0.100, STRAIN AMPLITUDE = .003 


3.0192 













0.50758 
0.80575 
Significance 
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TABLE 3.75: MEAN STRAIN = 0.100, STRAIN AMPLITUDE = .0025 


0.20743 
0.13162 | 0.13776 
0.87784 | 0.84225 
0 0 : 
> > 










059801 | 0.070682 
50.15 





4. Vote Tally 
Table 3.76 contains the vote count results for the 
fatigue life data. Since all the parameters were pooled, 
there is just one table. Data trends from the moment plots 
have eliminated Lognormal and Exponential; therefore, these 
distributions were eliminated in tallying the goodness of 


fit results. 


TABLE 3.76: FATIGUE LIFE VOTES 

[Moment Plots | 14 | 0 | 35. | 0 
[Goodness of fit | 26 | 0 | 37.) 0 
40] 0 | 72 | © 






Vote Totals: 


Based on the results of the random number testing of 
this method, a nearly 2:1 margin of Weibull votes over 
Normal votes is a strong indication that the samples come 


from a population with a Weibull distribution. 
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5. Weibull Zero-Shift 

Probability plots from the Weibull fit of each of 
the 16 sets of 20 samples are shown in Figures 3.77a through 
Be 77p. The plots are in the same order as the preceeding 
tables, so the mean strain and strain amplitude can be 
rererenced. When AGSS fits a sample to a Weibull 
distribution, it calculates the parameters of the Weibull 
distribution as well. To develop a model for the lives, 
these parameters are used. However, the fit for some of 
these samples can be improved by shifting the x-axis. 
Samples that have points dropping below the line at the 
lower tail indicate that a “zero shift" would improve their 


fet . 


WEIBULL PROBABILITY PLOT, N=20 WEIBULL PROBABILITY PLOT, N=20 





Figures 3.77a, b: Weibull Probability Plots 
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WEIBULL PROBABILITY PLOT, N=20 
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Figures 3.771, j, k, 1, m, n: Weibull Probability Plots 
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Figure 3.770, p: Weibull Probability Plots 


The plot of Figure 3.77g will be used as an example of how 
the zero-shift is accomplished. The curve of the left tail 
is extrapolated down to the x-axis as is done on Figure 
Caters The x-intercept is the amount that should be 
subtracted from the x-values, effectively shifting the zero 
to the intercept. Figures 3.78a, b, c, d, e, and £, show 
how the Weibull probability plots are affected by zero- 
shifts of -40,000, -42,000, -44,000, -46,000, -48,000, and - 


50,000 respectively. 
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Figure 3.78a, b, c, d, e, £: Weibull Probability Plots 


91 


Table 3.77 contains the goodness of fit statistics 


for the various zero-Shifts. 


TABLE 3.77: GOODNESS OF FIT VALUES FOR ZERO-SHIFTS 


Chi-Square | Kolm- Cramer-V.M. | Ander-Darl 
Smirn 


[None —~*4¢1 3487 
[=48,000 ~~ {0.52033_"[ 0.093617 [0.031944 0.27793 













From the table it appears that the shift that would produce 
the lowest values of the goodness of fit statistics is 
approximately -46,000. Figure 3.79 iS a graphical 


representation of Table 3.77. 


oz 


Chi-square 


Kolm.-Smirn. 


Cramer-V.M. 


Statistic Value 


~~~ Ander.-Darl. 
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Figure 3.79: Goodness of Fit Statistics vs. Shift Amount 


AGSS calculates shape and scale parameters of 2.2377 
and 16206 respectively when it fits the X7 data set to a 
Weibull distribution that is shifted to the left 46,000. To 
model the population of lives this data set is drawn from, 
these parameters would be used to generate lives from a 


Weibull distribution. 
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IV. SUMMARY AND CONCLUSIONS 
A. SUMMARY 

The moment plotting method on distributer 
characterization was developed by first deriving expressions 
for the coefficient of variation, measure of skewness, and 
measure of kurtosis, for the Normal, Lognormal, Weibull, and 
Exponential distributions. The derived expressions were 
plotted three ways as templates: (1) skewness’ versus 
coefficient of variation; (2) kurtosis versus skewness; and 
(3) kurtosis versus coefficient of variation. This produced 
three templates of moment curves. Before applying the 
method to fatigue data, the moment method was tested on 
random samples drawn from known distributions to determine 
whether or not it would correctly characterize the samples 
as being from the distribution they were drawn from. 

The first test used sample sizes of 1,000 and plotted 
sample moments on the templates. Points for the sample 
moments did not fall directly on any of the template curves, 
but trends were evident what were helpful in selecting the 
EGERPECE “AiStribucion. Template (1) exhibited horizontal 
trends in the data distribution for normal samples, while 
Lognormal and Weibull exhibited data trends with a 
Significant positive slope. Templates (2) and (3) were 
useful in distinguishing between Lognormal and Weibull 


distributions. Template (2) was distinguishing on the 
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negative skewness axis and template (3) possessed the 
largest spatial distance between Lognormal and Weibull. 

A vote counting method was also devised, where each 
plotted point voted for the distribution curve the point lay 
closest to on the template. Using this technique, the 
moment method correctly identified most of the _ sample 
distributions correctly; however, sometimes it did not and 
sometimes it was by a very narrow margin. When the sample 
Size was reduced to 100 and then to 20, the method grew more 
and more unreliable. A weakness of the moment plotting 
method was the fact that the moment functions forming the 
template were very near one another, and when the sample 
moments were plotted, they often spilled into the domain of 
more than one moment. 

The moment plots were useful, however, in eliminating 
distributions from consideration by observing the trends of 
the data and how they compared to the random samples 
initially tested. In most of the applications there was at 
least one distribution, and often times there were two, that 
could be ruled out as having trends unlike the data. This 
was the case in every one of the applications considered 
here. For instance, in the case of the helicopter loads 
data the Lognormal and Exponential distributions were 
dropped, reducing the task to choosing between the remaining 


two. 
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To improve the selection accuracy, goodness Olmems 
Statistics were coupled with the moment plotting method. 
Each sample was fit to a trial distribution using AGSS. For 
each fit, AGSS calculated the Chi-square, Kolmogorov- 
Smirnov, Cramer-von Mises, and Anderson-Darling goodness of 
fit statistics. The goodness of fit measures did not always 
agree among themselves, since they are measures of quite 
different things, and they were often times at odds with 
results from the locations on the moment plots, so the vote 
counting method was incorporated to make the final 
determination of the best distribution from the group 
remaining after the moment plotting method was used to make 
the first elimination(s). 

Applications were successfully made of the method to 
characterize the probability distribution functions of 
helicopter loads data, two different sets of tactical 
aircraft maneuver data, and a large set of experimental 
measurements of fatigue lives uSing uniaxial specimens. 

B. CONCLUSIONS 

Neither goodness of fit statistics nor the moment 
plotting method by themselves could consistently pick the 
correct distribution for a sample drawn from a known 
population. On the other hand when both were used together, 
employing the method outlined in the thesis, the success 


rate was 100%. 
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The method was applied to a broad range of fatigue data 
and found to be applicable. 

Applying the method to the several varied fatigue 
applications, there was considerable confidence in the 
Gistribution selected because of the preponderance of votes 
for that distribution by all the measures used in the final 
steps after the initial distributions were weeded out using 
the moment plotting method. 

This method, which employs both the trends that can be 
observed from the moment plotting method and the voting 
approach of the moments and the goodness of fit measures, 
seems to be considerably better than any of the measures 


used alone in a classical manner. 
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